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Abstract 



The classical and quantum dynamics of a particle trapped in a one- 
dimensional infinite square well with a time periodic pulsed field is inves- 
tigated. This is a two-parameter non-KAM generalization of the kicked rotor, 
which can be seen as the standard map of particles subjected to both smooth 
and hard potentials. The virtue of the generalization lies in the introduction 
of an extra parameter R which is the ratio of two length scales, namely the 
well width and the field wavelength. If R is a non-integer the dynamics is 
discontinuous and non-KAM. We have explored the role of R in controlling 
the localization properties of the eigenstates. In particular the connection 
between classical diffusion and localization is found to generalize reasonably 
well. In unbounded chaotic systems such as these, while the nearest neigh- 
bour spacing distribution of the eigenvalues is less sensitive to the nature of 
the classical dynamics, the distribution of participation ratios of the eigen- 
states proves to be a sensitive measure; in the chaotic regimes the latter being 
lognormal. We find that the tails of the well converged localized states are ex- 
ponentially localized despite the discontinuous dynamics while the bulk part 
shows fluctuations that tend to be closer to Random Matrix Theory predic- 
tions. Time evolving states show considerable R dependence and tuning R to 
enhance classical diffusion can lead to significantly larger quantum diffusion 
for the same field strengths, an effect that is potentially observable in present 
day experiments. 
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I. INTRODUCTION 



For several years now studies on quantized chaotic systems have increased significantly 
with the object of revealing quantum mechanical manifestations of classical chaos [p]||. The 
bulk of the work has used smooth Hamiltonian systems. If we start perturbing an integrable 
system, classical Hamiltonian chaos may develop through a gradual destruction of invariants. 
The celebrated Kolmogorov-Arnold-Moser (KAM) theorem gives conditions as to when a 
given tori would be only distorted. This scenario has been widely studied in two degree of 
freedom systems or two-dimensional area preserving maps. However, there are conditions 
upon which the KAM theorem rests that may not always be satisfied by certain systems 
of physical interest. In particular if the perturbation is not sufficiently smooth or even 
discontinuous the KAM scenario may break down. Large scale chaos may instantaneously 
develop in the system. One other way is that the KAM scenario fails when the unperturbed 
system is fully resonant, as in the Kepler problem. We deal in this paper with the former 
kind of non-KAM behaviour. 

We first discuss the prevalence of systems where such effects may be seen. The sim- 
plest systems where Hamiltonian chaos can develop is the so called 1.5 degree of freedom 
system, which are time-dependent one-degree of freedom systems. Thus consider the rotor 
Hamiltonian: 

H = f + f{t) V(9) 

where V(8) is an external potential that is periodic with period 2n, and f(t) is a periodic 
function of time with period T. A lengthening pendulum for instance may be the system 
under study. If V(8) is sufficiently smooth, the KAM theorem scenario combined with 
the Poincare-Birkhoff theorem provides the generic behaviour. The smoothness or at least 
continuity of V(6) is provided by the periodic boundary conditions in angular position of 
the rotor. Introducing discontinuous potentials will lead to delta function forces equivalent 
to walls of certain heights. 

This brings us to a natural class of systems where non-KAM behaviour will be the 
rule rather than the exception: externally forced particles in wells. This forms a broad 
class of systems which have evoked considerable interest and research since the development 
of quantum wells and dots. One of the experiments where quantum "scarring" of wave 
functions was reported involved resonant tunneling of a particle across a well in which there 
were external electromagnetic fields ||. 

In fact the simplest of such systems involve a particle in one-dimensional infinite square 
wells (1-d billiards) with time-dependent external fields. Consider as an example the Hamil- 
tonian: 

H = H + ecos(cjt)cos(27rx/A) (1) 

where H = p 2 /2 + V^a^a), describing such a particle. The potential V sq (x;a) is the 
confining infinite square well potential of width 2a, centered at the origin. Here e and A are 
field strength and wavelength of the external field which is being modulated in time with 
frequency uj. 
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It is easy to verify that the equations of motion are invariant under the following trans- 
formation: 



t — > u> t, p — > p/2au Q , x — > x/2a, 

and 

e — > e/(2au> ) , A — > A/2a, a; — > uj/ujq. 

Here the frequency c^o, which sets the new time scale, is arbitrary. Note that the new scaled 
variables and parameters are referred by the old symbols and they are dimensionless. Setting 
Uq = u in the above transformation, we have effectively two parameters : e and R = 2a/ X. 
Here R is the ratio of two length scales of the system i.e., well width and field wavelength. 
The presence of two competing length scales, provides a rich range of non-KAM behaviours. 
In particular if the dimensionless ratio R is a non-integer there is a possibility of observing 
non-KAM phenomena. 

Under the perturbation we can expect roughly that states whose absolute value of the 
initial momentum is less than y2~jej will be most affected. Thus low energy states will 
be most affected by the time-dependent forces. Fig. [I] shows the effect of the parameter 
R. While for R = 1 the system is essentially KAM and has KAM tori interspersed with 
resonances, any small deviation of R away from unity destroys low energy KAM curves and 
leads to increased chaos. Fig. |2| shows the fate of an individual KAM torus for which R = 1 
is a "bifurcation" point in parameter space and changes stability on either side. We expect 
such behaviour to be generic to a large class of similar systems and in this paper we will 
exhaustively study a "standard map" version of these systems. Just as the standard map 
provides an abstracted view of behaviour around nonlinear resonances we expect our model's 
analysis to provide such a view for these systems. 

Many models have been studied where the time dependence f(t) is a train of Dirac 
delta functions, the periodically kicked systems Q. The reasons are evident: they are the 
simplest Hamiltonian systems where many generic features of chaotic complex systems may 
be observed; and they allow a partial integration of the equations of motion, from kick to 
kick, enabling us to study iterative "mappings" rather than differential equations. Quantum 
mechanically the simplification enables us to partially integrate the Schrodinger equation 
and write the kick to kick propagator, or Floquet operator, analytically ||. 

An important paradigm in this class is the delta kicked rotor from which is derived the 
standard map ||[7|]. The dynamical richness of the classical system which obeys the KAM 
theorem and the consequent smooth transition to chaos is now well known, and fairly under- 
stood [jSHTOH . The corresponding quantum system has also been studied quite extensively as 



a model of "quantum chaos" (see |TTJ for an early review on this), and continues to provide 



an excellent model to numerically test our understanding of such systems ||12|| . 

The periodic input of energy into the system through the kicks can result in a diffu- 
sive increase of the momentum, and strong kicking strength can lead to unbounded energy 
growth classically. However, an important result on quantization is that the eigenfunctions 
(quasienergy states) are generically exponentially localized in momentum space, which sup- 
presses the momentum diffusion even in the highly chaotic regime. A plausible mechanism 
for this "dynamical localization" was suggested when an analogy was found El to Anderson 
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type localization of electrons in random on-site potentials []14j] . Experimental realizations of 



the delta kicked rotor, with cold atoms in pulsed standing laser fields [ 15fl , has confirmed 
the quantum suppression of diffusion. The localized states in unbounded momentum space 
results in quasi- independence of the quasienergies, and the random matrix properties |l6|,p| 
expected of quantized chaotic systems are not seen. For instance the nearest neighbour 
spacing distribution is Poisson rather than Wigner and the eigenfunction components are 
not gaussian distributed. It must be stated that most of these results are numerical and 
larger matrix calculations that are "more semiclassical" may show spectral transitions as 
the bulk of the eigenfunctions spread out and overlap with each other while the tails are 
still exponentially localized. Some evidence of this will also be present in this paper. 

Following our motivational discussion above, we replace the time dependence by a series 
of delta functions to facilitate the derivation of a map in which we can study non-KAM 



behaviour of the kind suggested above. Recently study of such systems has begun fll7| , |18 



For instance in |jnj it is shown that the quantum states are extended and delocalized in the 
highly chaotic (strong field) regime. In turn, the spacing distribution of the quasienergies, 
unlike the kicked rotor, follow the Wigner distribution. It is argued that the extended states 
do overlap and hence the corresponding quasienergies are not independent, resulting in level 
repulsion. Part of the present paper also critically examines these results. In [I8| a classical 
analysis of a generalized system has been carried out to understand the changes of stability 
that occur as a function of R and some of these results will be summarized below. 



II. CLASSICAL SYSTEM 

The system of interest is a particle inside the potential V sq (x; a) in the presence of a 
particular time-periodic impulse. We consider the Hamiltonian given by 

oo 

H = H + ecos(2rrx/\) ]T 8{n-t/T). (2) 

n=— oo 

Kick-to-kick dynamics of the particle immediately after each pulse can be described by an 
area preserving map which in dimensionless form is 

X n+1 = (-l) Bn {(X n + P n ) - Sgn(P n )B n } 
P n+ i = (-l) Bn P n + (K/2n)sm(27rRX n+1 ). 

Here B n = [Sgn(P n )(X„ + P n ) + 1/2] is the number of bounces of the particle between the 
walls during the interval between nth and (n + l)th kick, [...] stands for the integer part 
of the argument. The state of the particle just after the nth kick is now given in the new 
variables as X n ,P n . The sign of the momentum (±1) is given by Sgn(P n ). The following 
scaling relations are used to redefine the variables and parameters: 

X* P . = *f, K = ~~T~' (4) 
2a 2a aX A 

We note that \X n \ < 1/2 and effective parameters of the particle dynamics are K, the 
field strength and R, the ratio of the two length scales. The map in (JJ) shows the principal 
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features that we have discussed earlier in the introduction and may qualify as a "standard" 
map for such non-KAM systems. In fact there is a close relationship between this "well- 
map" (^) and the standard map itself that allows us to study it as a Generalized Standard 
Map (GSM) |H. This is the map: 

P n+1 = P n + (K/2tt) sm(2nRX n ) 

(5) 

X n+1 = X n + P n+ i (mod 1) 

which is defined on a cylinder (—00,00) x [—1/2,1/2). The well-map and the GSM differ 
only by boundary conditions i.e., the former and latter have reflective and periodic boundary 
conditions respectively. It is easy to see that this does not play any effective role, in the sense 
that the trajectories as evolved under the two systems can at most differ by a sign, depending 
on the number of bounces undergone. This fact simplifies considerably our analysis of the 
well-map. When R = 1, the GSM is the well studied and fairly understood standard map 
of the delta-kicked rotor. 

Dynamics of the GSM is highly chaotic and diffusive in the strong field regime (K ^> 1). 
In addition, it exhibits other interesting features like the development of chaos and hence 
diffusion even in the weak field regime (K < 1) when R 7^ j where j is a positive integer, 
see Fig. Note that the phase space portrait of both the well-map and the GSM are same. 
Such dynamical features are in fact common in non-KAM systems and one such situation 
is shown earlier. We can understand the development of chaos at low field strengths from 
the observation that the otherwise continuous map is discontinuous when R 7^ j. The KAM 
theorem does not hold for the discontinuous case and no smooth KAM tori exist in the phase 
space, however small K may be. In the absence of KAM tori, the phase space is chaotic and 
diffusive even in the weak field regime. More over, when R < 1/2 the GSM is a hyperbolic 



system. For a more detailed investigation of the GSM we refer the reader to |TS 



III. QUANTUM MECHANICS OF THE TRAPPED PARTICLE 

For the Hamiltonian which is periodic in time with period T(= 2tt/u), solutions of the 
Schrodinger equation satisfy the eigenvalue equation [|19| 

Utyi) = e-^ T /%> (6) 

where U is the one period time evolution operator. Here \tpj) and otj are the quasienergy 
states and quasienergies respectively. It is to be noted that the inner product of two ar- 
bitrary solutions of the Schrodinger equation for arbitrary time dependent Hamiltonian is 
independent of time due to hermiticity of the Hamiltonian. As a consequence of this, states 
correspond to and aj such that — aj 7^ nhu (where n is an integer) are orthogonal at 
a given time. If «j — aj = nhu, the states are degenerate and hence the quasienergies are 
uniquely defined with modulo Tiu. From the set of all orthogonal states we may write the 
general solution at a given time as — J2j Cjlrfj)- 
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A. Matrix Representation of U 



Periodically kicked systems are particularly easy to study since U can be written immedi- 
ately by integrating the Schrodinger equation between successive kicks. For the Hamiltonian 
H in Eq. (|2|) we have 

tf = exp{-*cos(^)}exp{-i^} (7) 

where k = eT/H. Note that the above time evolution operator is the quantum counter part 
of the well-map and not that of the GSM. The eigensystem in Eq. (|6]) may be solved by 
diagonalizing a matrix representation of U. The natural choice of basis for the [/-matrix is 
the eigenstates of the unperturbed Hamiltonian Hq: 

H \n) = E n \n) (8) 

where n = 1, 2, 3, The energy eigenfunctions and eigenvalues are 

J ^cos(ff), for nodd _ n 2 n 2 h 2 

W n > - j ^sin(ff), forneven ; En ~ 8a 2 " (!)) 

The unitary matrix U is then calculated as: 

U mn = (m\U\n) = (m\ exp{-iA; cos(27rx/A)}|n)e^ n2r = F mn e~ in2T (10) 

where we have defined an effective Planck constant as 

AT (id 



8a 2 

As the external field preserves parity we have 



where 



f 0, if m + n is odd 

I 27r \Q m - n ~ (— l)"Q m+n | , if m + 7i is even 



q 1 = r cos (w) e ~ ik ™(m d9 (13) 

J — 7T 



and 9 = nx/a. We note that Q\ is a Bessel function integral for integer R, while for non- 
integer R the integral constitutes a kind of "incomplete" Bessel function. Invoking the Bessel 
function J s (k) through the following identity 

oo 

e -ikcos0 = {-i) s J s {k) e~ isd 

s=— oo 

the integral can be evaluated as a series: 

oo 

Qi = 27rJ (k) 5 lfi + 2 £ {-i) s J s {k) C s (14) 
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where 



pit { (-l) l 2sRsm(sRw) r p / 1/1 

C s = / cos(Z0) cos{sR9) d9 = { (W^ 2 ' I0r sn T |£| . 
•>-■* I 7r, for sR = \l\ 

The relation J- S (k) = (—l) s J s (k) has been used in Eg. . Note that if R is an integer, 
Eq . (HH) simplifies to 

oo 

Q l = 271^ HY Js(k) S\ llsR (15) 

and a single term is picked out of the infinite series. 

The forms of Qi allow us to assess the fall of the matrix elements of the unitary matrix 
U. For integer R the unitary matrix can be essentially banded as the matrix elements fall off 
exponentially after a certain cut-off. For R = 1, as is well known and can be seen from above 
for I > k the matrix elements fall off exponentially, where / measures the distance from the 
diagonal. On the other hand when R is not an integer, apart from the Bessel function terms 
there are terms that are falling only polynomially in /. For instance when R = 1/2 we have 

Ql = m-i)' Mk) + (-i)<8 £ ( -' )a ;y 2) j a (fc). (i6) 

The infinite series gets effectively cut-off for s > k. The finite sum has terms that only 
decay as l~ 2 . Thus non-integer R values imply an important characteristic of the unitary 
quantum map: the polynomial fall of matrix elements, as opposed to the exponential fall 
characterizing integer R. In fact we may speculate whether non-KAM systems are always 
characterized by polynomially decaying matrix elements in the unperturbed basis. According 
to earlier studies eigenfunction localization crucially depends on the way in which matrix 
elements fall. 

When R = 1, the classical equivalence of the kicked rotor to the particle in a well was 
noted above. It is easy to see that quantum mechanically also the equivalence persists. The 
parity symmetry reduced rotor unitary matrix is identical to the well unitary operator in this 
case and hence odd states of the rotor correspond to the odd states of well, while the even 
states have a similar relationship. Thus all that is known for the quantum standard map, 
including exponential localization of eigenstates, may be carried over to the well system with 
R — 1. This allows us to address interesting questions of deviations from the standard map 
in a single model. 

The perturbing potential cos(27rx/A) preserves the parity of H , and hence U has the 
symmetry of parity. In what follows we consider only the states which have odd parity. In 
addition, the system has a spatial translational symmetry when R is an integer. Let us 
define a transformation for integer R as 

T/(X) = /((X + l/ J R)modl) (17) 

such that T R f(X) = f(X). T has the eigenvalues ft = exp(i^) where 1 = 0,1, 2...(R - 1). 
The commutation relation [U, T] = leading to T\ip) = For R = 2, (3i = ±1; in this 

case we consider only the states which correspond to (3i = 1. 
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Dimensionless quantum parameters k and r are related to the classical parameters 
through the relation K/R = 8kr. The semiclassical limit is k — * oo and r — > 0, such 
that fcr is fixed. Any arbitrary state of the system at a given time is = J2 n A n {t)\n) 

and its time evolution is given by A m (t + T) = Y. n U mn A n (t). 



B. Quantum Resonance 

Here we investigate if the parameter R has any effect on the important phenomenon 
of "quantum resonance". We notice that the unperturbed motion of the particle, given by 
the Hamiltonian H , between the kicks simply adds phase to the wave function components 
(when expressed in the unperturbed basis, as in Eq. fllPp). At resonance (r = 2%), the 
unperturbed motion between the kicks is absent. In this case, without loss of generality, the 
time evolution of an arbitrary state of the system is 

|#(t)) = e -^cos(2WA)i|^( )) (18) 

and thus |\I/(t)| 2 = |\I/(0)| 2 . Note that here t is the number of kicks. The kinetic energy of 
the particle is then 

f |^(0)| 2 sin 2 ^)dx|. (19) 



'2ttH\ 2 r 



V A 

In the limit t — > oo energy grows quadratically with the number of kicks. If |^(0)) = |n), 
i.e., the initial state is one of the unperturbed state itself, the energy is purely quadratic. In 
fact, the energy can be found exactly as 

hV) = hiO){l+(—)\2-A)} 120) 



n 



where 



.4 



sin(27rfl) 
I 

\n+l 



ifn = 2R 



Since A ^ 2, we observe that the quadratic energy growth is unaffected by the length 
scale ratio R. Numerically we have found that this behaviour is seen when r is rational 
multiples of 2n also. Thus the quantum resonance phenomena of the well system is very 
similar to that of the kicked rotor PJ2"U|. It is to be noted that resonance is a non-generic 
pure quantum phenomena and no correspondence to it can be seen in the classical system. 
In the context of a particle in a well, quantum resonance may lead to enhanced ionization 
in a finite well. 
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IV. RESULTS 



Having given sufficient description of the system under investigation, here we analyze 
quasienergy states and quasienergies of the generic quantum system (r is irrational multiples 
of 2-71") in relevant classical regimes. On taking a truncated A dimensional Hilbert space 
spanned by the first iV unperturbed basis states which belong to odd parity, diagonalization 
of the matrix U mn gives the eigenstates such that = Ysni>n\ n )- We consider 

only states that are "converged" in the sense that they are independent of the truncation 
size A. Thus the states we are interested in belong to that of the infinite Hilbert space; 
they are states of the infinite cylinder and not of a truncated cylinder, or torus. The 
last distinction becomes important as quantum states that belong to the cylinder can have 
completely different localization features than those that belong to a truncated cylinder. As 
we are interested in a particle in an infinite potential well, such a truncation lacks physical 
meaning. 



A. Localization measures of eigenstates 

Localization can be measured using a unified quantity, the Renyi Participation Ratio £ q 

/ \ 1/(9-1) 

^=(X>n| 2 j (21) 

of which the entropy and participation ratio (PR) are special cases. In our analysis we 
first use a normalized information entropy as a measure of localization of states, and this is 
defined as: 

^E^jhW'M*.!' (22) 

It is easy to see that S = ln£i/ In (A/2). This measure compares the entropy to that of the 
eigenfunctions of A x A matrices belonging to the Gaussian Orthogonal Ensemble (GOE) 
which is approximately ln(A/2). The GOE is relevant to time reversal symmetric systems 
such as we are considering. 

First we calculate a gross measure of localization in a given spectrum by averaging over 
all converged states. We set criteria for the states to be converged so that the states belong 
to the cylinder, or are at least very close to states that belong to the cylinder. In all the 
following cases, the eigenvalues are converged in modulus to unity to within 0.0001 or better. 
Fig. |] shows the average entropy as a function of R. For small K (< 1), the oscillations 
are qualitatively similar with distinct entropy minima at integer R and maxima at around 
half-integer R. This may provide a simple mechanism to experimentally control the extent 
of localization. The information entropy is of course basis dependent; the unperturbed basis 
we use is a useful one as it has information about localization in the momentum. 

Naturally, the minima in entropy is expected to have strong associations with the pres- 
ence of stable regions in the classical phase space. Of special significance are KAM tori in 
phase space, as these structures are complete barriers to classical diffusion in momentum. In 
spite of the fact that in the classical system all the KAM tori break up in the standard map 
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(i? = 1) at K = 1, we observe a minimum entropy. This is due to the presence of cantori 
which are partial barriers for chaotic orbits and suppresses global diffusion. For non-integer 
R a complex phase space picture emerges and has been discussed in |18| . Maximum entropy 



around half-integer R is the classical parametric regime where the discontinuity is maximum 
corresponding to the maximum chaos assisted diffusion. 

For large K (= 10), oscillations in entropy are still present while there is apparently 
complete chaos for all relevant R values. We can understand these oscillations as due to 
the strong correlation between the localization of eigenstates and classical diffusion coeffi- 
cient. For R < 1/2 the semiclassical parameter k = K/(8Rt) is large, yet there is increased 
localization of states due to limited classical diffusion, presumably due to the presence of 
cantori. For the kicked rotor the exponential localization length was found to be propor- 
tional to the classical diffusion coefficient This was found by numerical experiments 
and is supported by certain qualitative arguments. We are now in a position to examine 
the relationship between quantum localization and classical diffusion in the context of the 
particle in a well, wherein we have the freedom of another control parameter, namely R, 
with which to vary the classical diffusion. 

Instead of studying localization lengths we study here the measures of localization such 
as the entropy or the PR. We study the PR more closely rather than the entropy. In chaotic 
regimes we have numerically ascertained that exponential of the entropy is proportional 
to the PR, as shown in Fig. || The relationship between the localization length hitherto 
calculated for the kicked rotor and the PR calculations we present will need more detailed 
study, but we expect them to be roughly proportional to each other. In fact if we assume a 
fully exponentially localized state with \if) n \ ~ exp(— \n — n |//oo), then the PR is 

& = (EWT 1 = 2 ^oo. (23) 

n 

We recapitulate the argument connecting classical diffusion and the localization length 
for the specific system we are considering, as there are difference in factors. On considering 
time evolution of an initial state, kinetic energy diffuses for a certain time t c and then attains 
quasiperiodic saturation. The number, n c , of unperturbed states that are excited during the 
time evolution is related to the critical time by the diffusion equation: 

7T 2 h 2 n 2 c = D CL t c ; (( Pt - Pof) = D CL t (24) 

where D Ch is the classical diffusion coefficient in momentum and (..) represents the ensemble 
average. Here the momenta and the diffusion coefficient have dimensions and we have 
taken a = 1/2. The critical time being the Heisenberg time relevant for n c equally spaced 
eigenstates, we get t c ~ n c T '/2tt. If the average localization length (l^) is also n c , we obtain 
the relation : 

(£ 1 ) = 2(l OB ) = ^D{K,R) (25) 

where r is the dimensionless effective Planck constant defined in Eq. ([11]) and a is a constant 
whose value has been numerically determined as 1/2 for the standard map [2^]. D(K, R) is 



the dimensionless diffusion coefficient which one will get from using the dimensionless maps 
Eq. (0) or Eq. (R). The dependence on both K and R is emphasized. 
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In Fig. |6] we show the average PR and the scaled diffusion coefficient according to the 
relation Eq . (|25|) . We see that the relation derived above holds in some parameter regions 
while it picks up only qualitative features of the oscillations in the others. In particular 
the relation seems to hold for R < 1/2 when the classical system is hyperbolic as well 
as around R — 1. The deviations from the relation fl25|) might be due to fluctuations of 
the state components in the unperturbed basis (one such case is shown in Fig. |TT]). These 



fluctuations may lead to different scaling behaviour between the average PR and the classical 
diffusion coefficient. However more detailed investigations are needed to make any positive 
statements. The sharp deviation for R = 2 can be accounted for as due to the presence of 
an extra quantum symmetry discussed above. 

Following our study on average PR and its scaling with the classical diffusion coefficient, 
we may then enquire about how the PR itself is distributed in a given spectrum if the 
average reflects the general behaviour. We find that when the classical system is chaotic, the 
distribution of a normalized quantity y = ln£f / (ln^ 1 ) (this is similar to the distribution 
of the entropy due to the linear relationship exhibited above) is nearly normal as seen in 
Fig. |7|. This may be attributed to a realization of the Central Limit Theorem. However the 
PRs and Inverse Participation Ratios (IPRs) themselves are not normal. Their distributions 
may be got by assuming that the distribution of y is normal. Thus the PRs are distributed 
according to the lognormal distribution [p3|: 



^- ^■■W H-^tlf?- 1 " (26) 




where a 2 is the variance of y. As an immediate consequence, distribution of IPRs is also 
lognormal. Distribution of such localization measures are much of significance. Recently the 
distribution of IPRs have been exploited to show that the distribution of resonance widths 
in wave-chaotic dielectric cavities is lognormal Jj4[ . 



When K is small (< 1), the classical motion is nearly regular for R — 1, while chaotic 
for R < 0.5. However the time scales for classical diffusion is large making the observation 
of R effects on quantum dynamics hard to discern. For instance the nearest neighbour 
spacing distribution may remain very close to the Poisson distribution. In such a situation 
we find that the distribution of the PRs provides a positive litmus test. In Fig. [| such an 
example is shown, wherein even for small field strengths the effect of R is clearly visible as 
a tendency for y to be normally distributed. This is an indication of the "derealization" 
that is taking place in the eigenfunctions. This derealization is limited in the sense that 
while the eigenfunctions remain square integrable there is more spreading out in the bulk 
part of the states. Thus we may conclude that distribution of the localization measures is a 
sensitive quantity in chaotic quantum systems. 

Time evolution of non-stationary states must reflect the properties of the stationary 
states and is also of importance in the context of experiments. Here we have studied the 
diffusion in kinetic energy of a state \^>) which is initially the ground state of the unperturbed 
system. We illustrate with one example wherein for a fixed classical parameter K, the effects 
of non-integer R is seen clearly for a given r value. Thus tuning R essentially tunes A since 
a is fixed through the relation (|TTp. In Fig. |D| scaled kinetic energy (P 2 ) = (^\P 2 \ty) is 
shown as a function of time (number of kicks) for a small value of K corresponding to a 
small classical field strength e. We note that while the quantum diffusion saturates at a 
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much higher value for R = 1.5, compared to R = 1, the actual classical field strength e 
(from Eq. (^)) is smaller by a factor of 1.5. For comparison is shown another integer case, 
R = 2, wherein the classical diffusion is smaller compared to R = 1. 



B. Eigenvalues and Eigenstates 

It is clear from our earlier observations that the states are more localized in the regular 
or mixed regimes of the classical system while less localized (or delocalized) in the chaotic 
regimes. The degree of localization is also controlled by the ratio of the length scales and 
complexity of the classical phase space is reflected in the localization measures. Here we 
look at the quasienergies and the corresponding states more closely. 

In Fig. |1(] we have shown the nearest neighbour spacing distribution of the quasienergies 
for various parameters. The first row and the last column of the catalogue correspond to 
classically chaotic regimes and the rest belong to the regular/mixed phase space regimes. 
In regular/mixed regimes where the states are highly localized, the spacing shows excellent 
agreement with the Poisson distribution. On the other hand, in chaotic regimes the spacing 
agrees well with the Poisson distribution except at small spacings. This is due to the fact 
that bulk part of the eigenstates are delocalized and they overlap each other. However, the 
tail part of the states are exponentially localized and the degree of overlap is not significant 
enough. We also notice that the spacing distribution is only slightly sensitive to the nature 
of the classical dynamics in the case of the unbounded kicked rotor or the well, at least in 
the parameter regimes we have investigated. In such situations, as we have demonstrated 
earlier, the distribution of PRs is a good measure to distinguish the chaotic quantum systems 
from the regular systems. 

Our extensive calculations of the eigenstates in chaotic regimes show that, in general, it 
is hard to qualitatively differentiate the states corresponding to non-integer R values from 
the rotor (R = 1) states as far as their localization behaviour is concerned. In particular it 
is not easy to distinguish the emergence of non-exponential tails unequivocally. However, 
we found that eigenstates corresponding to non-integer R values generally have more fluctu- 



ations compared to the rotor states; this is illustrated with some examples in Fig. 0]. The 
fluctuations are closer to the RMT predictions in the case of non-integer R values and is 
shown further below. 

Recently there have been studies on the special case (R = 0.5) of the system (|2|), with 
the motivation of revealing quantal behaviour of non-KAM systems |]17|. It was observed 



that the quasienergy states are "extended" in the unperturbed basis and as a result the 
spacing distribution was shown to be Wigner distributed. At this juncture we would like to 



compare our results with certain aspects of this work. In JT7[], the eigenstate shown in the 
highly chaotic regime (K = 50, N = 1024; we have not been able to ascertain the value of 
r used in this work) does not appear to belong to the unbounded phase space as it spreads 
all over the basis. Thus while states such as these may belong to some truncated dynam- 
ical system, it does not belong to the infinite Hilbert space of the well system. Increasing 
the dimensionality of the matrix used will modify such states; in short they are not con- 
verged. As we demonstrate below, unconverged or poorly converged states may mislead us 
in understanding the spectrum. 
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Large K implies large k for given R and r, and hence our calculation demands bigger 
dimensionality iV of the truncated Hilbert space, since the PR is roughly increasing as k 2 . 
Although we take N = 2000, getting a good number of converged states is problematic. 
We pursue the spacing distribution with different convergence criteria for the numerically 
obtained states. The convergence criteria uses partial sum of the state components : 



M 

|2 



{Sum} M = ]T |V„| 2 ! with M<N. (27) 



n=l 



For a well converged state we expect that {Sum} M m 1, even for M <^ N . We denote by 
Nm the number of converged eigenstates whose {Sum} M is greater than Sm (an arbitrary 
number close to, but less than, unity) for a fixed value of M. Thus the convergence criteria 
is characterized by M and Sm- 



In Fig. [12|, we have shown the spacing distributions with different criteria for two cases. 
In both the cases transition to Wigner distribution is evident as the convergence criteria is 
relaxed. The unconverged or poorly converged states do not belong to the physical system 
of our interest and the corresponding quasienergies follow the RMT prediction. Obviously, 
reliability of the result is more in the top plots where the spacing shows neither Poisson nor 
Wigner distributions. Though the tail part shows the Poisson behaviour there are significant 
discrepancies in the small spacing. A more correct picture may be closer to the scenario of 



the chaotic regimes presented in Fig. 10 



Shown in Fig. [13] are few "well converged" states, with more stringent convergence 
criteria (M = 1600, Sm — 0.9999). With this criteria we have only Nm = 12 and 4 for 
R = 1 and 0.5 respectively. The state components exhibit strong fluctuations in the basis. 
Here again it is hard to differentiate the two cases qualitatively. The states corresponding 
to R = 0.5 also appear to have exponential tails. To see the distribution of the state 
components, we introduce a variable r] n = l^nlVl^nl 2 where over bar stands for average over 



the state components such that r) — 1. As seen from Fig. [L4|, the cumulative distribution 



of rj for both the cases have very similar behaviour. Considerable deviations from the RMT 



predicted cumulative Porter-Thomas distribution, I(rj) = Erf (y r//2j, may be attributed to 
the localization of the states. However, the distribution corresponds to R = 0.5 tends to be 
closer to the RMT predicted behaviour. 



V. SUMMARY AND CONCLUSION 

A particle inside a one-dimensional infinite square well potential in the presence of a time 
periodic pulsed field is examined both classically and quantum mechanically. This simple 
model can be seen as one generalization of the kicked rotor, or the standard map. A variety 
of classical dynamical features emerge from the nature of the ratio, R, of the two competing 
length scales (the well width and field wavelength). Many of the dynamical features so 
observed are generic to a wide class of systems of substantial current interest viz. externally 
forced particles in wells. 

It is shown that when the length scales do not match, even in perturbative regimes dy- 
namics can be increasingly complex wherein all the KAM tori in phase space break up. As a 
result the transition to chaos is abrupt, a typical scenario of non-KAM behaviour. Quantum 
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mechanically imprints of such a transition is seen as a spread in bulk part (delocalization) 
of the eigenstates. Thus we realize the length scale ratio R as a control parameter for the 
localization in the weak field regime. 

On increasing the field strength, chaos assisted diffusion takes place in momentum. From 
earlier studies on the kicked rotor it is known that the average localization length of the 
eigenstates is directly proportional to the classical diffusion coefficient. We have shown that 
in our generalization of the kicked rotor also, this result grossly explains the localization 
behaviour of the eigenstates through the classical transport property. Thus the kicked rotor 
continues to serve as a useful model in understanding physical phenomena exhibited by a 
larger class of systems. 

We have observed, as expected, that in the regular case the nearest neighbour spacing 
distribution of the quasienergies show good agreement with Poisson distribution. We have 
presented evidence to support that, in highly chaotic regimes spacings show some deviations 
from the Poisson distribution though the corresponding eigenstates belong to an unbounded 
phase space. Limited overlap of the eigenstates results in such deviations. However, spacing 
does not show the RMT predicted Wigner distribution as was claimed in an earlier study 



Yt\ . The earlier result is attributed to lack of converged states which make up the statistics. 

While spacing is not much sensitive to the classical chaos, the distribution of participation 
ratios of the eigenstates is shown to be a good measure to distinguish chaotic quantum 
systems from the regular ones. Quantum mechanically, chaotic regimes are characterized 
by a lognormal distribution of the participation ratios. In addition to the above generic 
quantum features, we have also studied non-generic phenomena like "quantum resonance". 
In resonance condition, kinetic energy of the particle grows quadratically with the number 
of kicks. This unbounded energy growth is not affected by the length scale ratio and can 
enhance ionization in the finite well system. 

As far as experimental realizations of this work is concerned, perhaps both quantum 
wells H and cold atom experiments are possible candidates. As suggested above, the 
R effects may be best observed at small field strengths and for R > 1. Further work is 
underway exploring the nature of localization in such systems, including a bounded version 
of the generalized standard map. 
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FIGURES 




and to = 1. The lower momentum region is increasingly chaotic when the length scales do not 
match. 
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FIG. 2. Shown are orbits of Fig. [I] having identical inital conditions. The initial condition 
corresponds to a KAM torus in the lower momentum region for R = 1 (the negative momentum 
region is not shown here). Note the abrupt change in the stability and the non-generic features of 
the resultant phase space structures. 
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R = 0,95 R = 1 R = L05 




FIG. 3. Phase space portrait of the GSM with K = 0.3. For R = 1, the dynamics is nearly 
regular wherein many smooth KAM tori are seen. When R departs from unity the dynamics is 
increasingly complex and no KAM tori are seen. This may be compared to the lower momentum 
region in Fig. |]. 
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FIG. 6. The average PR (•) and the scaled classical diffusion coefficient (o) are plotted as a 
function of R for the case K = 10, r = 0.1. The dotted line is the scaled coefficient calculated using 
up to the second order time correlation. Higher order time correlations are insignificant since the 
classical system is highly chaotic. 
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FIG. 7. Probability distribution of y, the normalized log of the PR, for the kicked rotor 
(R = 1) in the chaotic regime. Here we have taken K = 10 and r = 0.025 (x), r = 0.05 
Smooth curves are corresponding gaussian distributions. 
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FIG. 8. Probability distribution of y, the normalized log of the PR, for the case 
K = 0.1, r = 0.001 (first row) and for K = l,r = 0.01 (second row). Smooth curves are 
corresponding gaussian distributions. Note the sensitivity of these distributions to the classical 
dynamics. 
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FIG. 9. Shown is the scaled kinetic energy (P 2 ) of a state, which is initially the ground state 
of the unperturbed system, as a function of time. Here the parameters are K = l,r = 0.01 ; 
R = 1 (solid line), R = 1.5 (dots) and R = 2 (dotted line). Effect of non-integer R is remarkably 
seen in the evolution as the kinetic energy of the quantum particle saturates at much higher value 
compared to integer R cases. 
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FIG. 10. Nearest neighbour spacing distributions of 1000 quasienergies for (a) 
K = 0.1, r = 0.001; (b) K = l,r = 0.01; (c) K = 10, r = 0.1 with R = 0.5,1,1.5 (top to 
bottom) and N = 1200. Smooth curves are Poisson distributions. Note the relative insensitivity 
of these distributions to the classical dynamics. 
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FIG. 12. The nearest neighbour spacing distributions for the case K = 50, r = 0.1. Smooth 
curves are the Poisson and Wigner distributions. Convergence criteria is relaxed as we move from 
top to bottom. A "spectral transition" is observed. 
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FIG. 14. Collective cumulative distribution of the components of the states shown in Fig. 13. 
Dotted curve corresponds to R = 1 while dashed curve corresponds to R = 0.5. Solid curve is the 
cumulative Porter-Thomas distribution. 
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